Asymptotic expansions for high-contrast elliptic 

equations 
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1 Abstract 

In this paper, we present a high-order expansion for elhptic equations in high- 
contrast media. The background conductivity is taken to be one and we assume 
the medium contains high (or low) conductivity inclusions. Wc derive an asymp- 
totic expansion with respect to the contrast and provide a procedure to compute 
the terms in the expansion. The computation of the expansion does not depend 
on the contrast which is important for simulations. The latter allows avoid- 
ing increased mesh resolution around high conductivity features. This work is 
partly motivated by our earlier work in [22| where wc design efficient numeri- 
cal procedures for solving high-contrast problems. These multiscale approaches 
require local solutions and our proposed high-order expansion can be used to 
approximate these local solutions inexpensively. In the case of a large-number 
of inclusions, the proposed analysis can help to design localization techniques 
for computing the terms in the expansion. In the paper, we present a rigorous 
analysis of the proposed high-order expansion and estimate the remainder of it. 
We consider both high and low conductivity inclusions. 

2 Introduction 

The mathematical analysis and numerical analysis of partial differential equa- 
tions in high-contrast and multiscale media are important for many practical 
applications. For instance, in porous media applications, the permeability of 
subsurface regions is described as a quantity with high-contrast and multiscale 
features. A main goal is to understand the effects and complexity related to 
this multiscale variation and high-contrast in the coefficients. This is specially 
important for the computation of numerical solutions and quantities of interest. 
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Many tools and methods have been developed and used to study high contrast 
problems. We mention few recent works where numerical methods are designed 
that target problems with difficult variation in the coefficients. The numeri- 
cal analysis of these methods require accurate descriptions of the variation of 
the coefficients. In 0, 0, H E [ll [3 E IH S, |33| multiscale methods for 
problems with high-contrast coefficients are described. For domain decompo- 
sition methods with discontinuous coefficients we mention 15. 28l 30l 34. SSj 
and references therein. Some domain decomposition techniques for multiscale 
partial differential equations with complicated variations in the coefficients are 



developed in [16|, UJ, |22|-|2J, 12; 



UlllZllaal among 



33j . Additionally, we mention 
others, which focus on multilevel methods targeting problems with high-contrast 
and discontinuous coefficients. Numerical methods based on different asymp- 
totic analysis are detailed in 0, H [H , 21, 27 1 . 

In the case of elliptic problems with oscillating coefficients and bounded con- 
trast, multiscale methods and homogenizations techniques have been success- 
fully applied to study solutions of elliptic differential equations (c.f., [s, 2^ 3^ 
and references therein). Many homogenization and multiscale methods start 
with the derivation of an asymptotic expansion for the solution of the partial 
differential equation. The expansion is written in terms of the problem pa- 
ramcter(s): e.g., the period in the case of oscillating periodic coefficients. The 
asymptotic expansion is then used to study the problem at hand. 

In this paper and in the same spirit, we derive asymptotic expansions for 
the solutions of elliptic problems with high-contrast. In this case the parameter 
to be consider is the contrast in the coefficient. In particular, we consider the 
problem 

- div(K(a;)Vu) = /, in D (1) 

with Dirichlet data given by it = g on dD. We assume that k{x) > 0. The con- 
trast in the coefficient, ||k||^oc(£))||k^^||^oc(^), is the important parameter con- 
sidered here. For the analysis, we consider a binary media k{x) with background 
one and with multiple (connected) inclusions. We consider high-conductivity 
inclusions (with conductivity rj) and low-conductivity inclusions (with conduc- 
tivity 1 /r]) . Wc derive expansions of the form 



= 77W_l + Uq + 



1 



(2) 

and © 
(3) 



In the case with only high-conductivity inclusions we have w_i 
reduces to 

1 1 

u,, = ua H — ui H — + 

rj 

In the presence of low-conductivity inclusions we may have U-i ^ depending 
on the support of the forcing term. We mention here that, in order to derive 
the expansions, we use the weak formulation associated to ([IJ. In this case, 
using the integral formulation has the advantage that the boundary, interface 
and transmission conditions are self-revealing. 
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The asymptotic problems for the case of only high-conductivity or only low- 
conductivity are studied in detail. For the study of the high-conductivity in- 
clusions asymptotic problem we use harmonic characteristic functions. These 
functions are defined as being constant inside the inclusions and harmonic in 
the background domain. The asymptotic solution can be obtained by solving a 
Dirichlct problem in the background domain and a finite dimensional problem 
in the space spanned by the harmonic characteristic functions. The solution of 
the finite dimensional problem gives a closed formula for the constant values 
of the limit solution inside the high-conductivity inclusions. The resulting sys- 
tem can be large, in general, and one can consider some localization techniques 
(c.f., ii). 

The asymptotic problem and approximations of the asymptotic problem, in 
the case of high-conductivity inclusions, have been studied in the literature. 
We mention 0, 0, where a discrete network approximation is considered for 
the problem of computing the effective conductivity of high-contrast, randomly, 
and densely packed composites with high-conductivity inclusions. The network 
approximation depends on the geometry of the inclusions and the behavior of 
the solution between nearby inclusions. The authors can localize the interaction 
of high-conductivity inclusions using graph-theoretical concepts. Furthermore, 
they propose a finite element approach for solving the resulting system and iden- 
tifying the first order approximation of the solution. In Q the authors study 
the homogenization of the asymptotic problem in terms of geometric parameters 
such as the shapes of inclusions and the distance between the inclusions. They 
develop an asymptotic analysis for periodic structures with absolutely conduc- 
tive square inclusions. The small scales considered here are the period of the 
structure and the distance between inclusions. We refer to the works and 
references therein. 

We also write a low-conductivity asymptotic problem valid only in the case 
where the forcing term vanishes inside the low-conductivity inclusions. This 
problem is important in flow applications where low conductivity regions repre- 
sent shale regions and can substantially alter the overall flow behavior. To our 
best knowledge, this problem is not extensively studied in the literature. The 
asymptotic solution can be obtained by solving: 1) a mixed boundary condition 
problem in the background domain, and then, 2) a Dirichlet problem inside the 
inclusion with zero forcing term and the Dirichlet data from 1). 

We show how to obtain all the coefficients in the expansions. The procedure 
to compute the coefficients, coincide with a Dirichlet-to-Neumann procedure 



as in Domain Decompositions Methods; see [28|, |35|. For the case of high- 
conductivity inclusions, the Neumann problems are solved in the interior inclu- 
sions and therefore, a compatibility condition needs to be verified. Obtaining 
the right balance of fluxes for the compatibility condition involves the solution 
of a finite dimensional problem in the space spanned by the harmonic charac- 
teristic functions mentioned above. For the case of low-conductivity inclusions 
the Neumann problems are solved in the background domain and no compati- 
bility condition in required. The expansions derived in this paper are proven to 
converge in H^{D) for high-contrast bigger that a certain constant. This con- 
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stant depends on the domains representing the inclusions and the background 
domain. Asymptotic expansions in the presence of boundary intersecting in- 
clusions can be derived and analyzed using similar arguments. The presence 
of low- and high-conductivity inclusions can be also analyzed. More general 
geometrical configurations and partial differential equations can be studied as 
well. 

Having a practical procedure to compute the next leading order terms in ([2]) 
is useful for applications. For instance, the quantity ui in ([3]) may have con- 
siderable contribution to some quantity of interest in some regions, e.g., in the 
velocity k|Vm^| expansion, the solution is multiplied by ry. A high-order expan- 
sion is also useful when constructing multiscale and multilevel methods. The 
expansion ([2]) can be used to construct multiscale finite element basis functions; 



see 



12l . |20| . Such an expansion will allow the construction of basis functions 
independent of the contrast and depending only on the limiting problem. In the 
case of expansion ([2]), the asymptotic problem depends only on the geometry 
configuration describing the inclusions (see Section l4.3.1|) . These basis functions 
will capture the effect on the solution of the geometric arrangement describing 
the conductivity. The next order terms in ([2]) can be used to construct correction 
terms to account for the effect of the contrast in the coefficient. Expansion ([2]) 
can be used to improve existing advanced multiscale finite element techniques 
for a better sub-grid capturing; see l3, l^l . Fast numerical upscaling techniques 



can also be constructed with the first order terms of © or ([3]). See (2l| where 



the authors develop fast numerical upscaling methods based on some asymp- 
totic analysis. We also mention [3, H^l where numerical approximations are 
designed using asymptotic analysis. 

The rest of the paper is organized as follows. In Section [3] we recall the 
weak formulation of (H]). In Section!?] we derive the expansion for the case of 
high-conductivity inclusions. Wc study the asymptotic problem and the conver- 
gence of the expansions. Section [5] is dedicated to the case of low-conductivity 
inclusions. In Section [6] we consider the case with low- and high-conductivity 
inclusions and in Section [7] we make some conclusion and final comments. 



3 Problem Setting 

Let D CM!^ polygonal domain or a domain with smooth boundary. We consider 
the following weak formulation of Find u E H^{D) such that 



a{u,v) = f{v) yv^HliD), 
u = g on dD. 



(4) 



Here the bilinear form a and the linear functional / are defined by 

a{u,v)= / k(x)Vu(x) ■ Vi)(x) yu,v e H^{D) (5) 

Jd 

f{v)= I f{x)v{x) WveH'oiD). (6) 
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We assume that D is the disjoint union of a background domain and inclusions, 
D = DqU {U^^iDjn)- We assume that Dq, Di, . . . , Dm, are polygonal domains 
(or domains with smooth boundaries). We also assume that each Dm is a con- 
nected domain, m = 0, 1, . . . , M. Let Dq represent the background domain and 
the subdomains {Dm}m=i represent the inclusions. For simplicity of the pre- 
sentation we consider only interior inclusions. See Figure [T] for two dimensional 
illustrations. 




Figure 1: Examples of geometry configurations with interior inclusions. 

Given w G H^{D) we will use the notation w^^^\ for the restriction of w to 
the domain I?„i, that is, 

lyM^M;!^,^, m = 0,l,...,M. 

4 Expansions for high-conductivity inclusions 

In this section we derive and analyze expansions for the case of high conductivity 
inclusions. For the sake of readability and presentation, we consider first the case 
of only one high-conductivity inclusion in Scction [4.1l and study the convergence 
of this expansion in Section 14.21 We present the multiple high-conductivity 
inclusions case in Section H31 where we describe the expansion and analyze its 
convergence, following the structure presented in Sections 14.11 and 14.21 

4.1 Derivation for one high-conductivity inclusion 

Let K be defined by 

'^^"^^{t .eD,^=D\D^, 

and denote by the solution of the weak formulation ([4]). We assume that 
Di is compactly included in D [Di C D). Since Urj is solution of (j4]) with the 
coefficient (O, we have 

/ \/u,,-Vv + r] I \/u,,-Vv^ I fv \/veH^{D). (8) 

JDa J Di Jd 
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We seek to determine {uJ^q ^^i^) such that, 

11 °° . 

M„ = Mo + -"1 + -^U2 H = 77"*Mi, (9) 

' ' 1=0 

and such that they satisfy the foUowing Dirichlet boundary conditions, 

uq = g on dD and = on dD for i > 1. (10) 
We substitute ^ into ([S]) to obtain that for all v G Hq{D) we have, 

OO n OO n n 



or 



Vt; = / /«. (11) 



7? / Vuo • Vw + V 77"M / Vuj ■Vv+ Vm,+i ■ Vu ) = / 

Now we collect terms with equal powers of rj and analyze the resulting subdom- 
main equations. 

4.1.1 Term corresponding to 7/ = rj^ 

In (|lip there is one term corresponding to rj to the power 1, thus we obtain the 
following equation 



/ Vuq • Vw = for all v e (D). 

JDi 



(12) 



In the general case, the meaning of this equation depends on the relative position 
of the inclusion Di with respect to the boundary. It may need to take the 
boundary data into account. Since we are assuming that Di C D, we conclude 
that VitQ^"* = in Z?i and then u\f' (the restriction of uo on Di) is a constant. 

4.1.2 Terms corresponding to 77*^ = 1 

Equation pT|) contains three terms corresponding to rj to the power 0, which 
are: 

/ VTio • Vu + / Viii • Vw = / fv for all v e HI{D). (13) 

J Do J Di JD 

Let 

Vconst = {w e Hq{D), such that w'^' = is constant}. 

If we consider z S Vconst in equation we conclude that uq satisfies the 
following problem, 

/ VU0-VZ=) / V7to-Vz= / fz Vconst (M) 

/ Jdo Jd 

uo = g on (15) 
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The problem is elliptic and it has a unique solution. To analyze this 
problem further, it is natural to define a harmonic characteristic function xdi S 
H^iD) such that 

xL'! = l ini^i, 
and the harmonic extension of its boundary data in Dq is given by 

J Do 

X^d\ = 1 on dD^, (16) 

=0 on aL*. 

To obtain an explicit formula for uq we will use the following facts: {i) prob- 
lem (jl4p is elliptic and has a unique solution, and (m) a property of the harmonic 
characteristic functions described in the Remark below. 

Remark 1 Let w be a harmonic extension to Dq of its Neumann data on DDq . 
That is, w satisfy the following problem, 

/ Vw-Vz^ Vwnoz for all z € (Do). 

Jdo JdDo 

Since xDi ~ on dD and xDi = 1 on dDi , we readily have that 



^XDi ■ Vw = / Vw ■ nxDi = ( / ' ) ^ M / ' 

JdDo \JdD / \JdDi 

and we conclude that for every harmonic function on Dq, 



VxDi ■Vw= Vw no. (17) 

Ai JdDi 

In particular, taking w = xDi we have: 

I\'^XdA^=I Wxd,?=I VxD,-no. (18) 

Jd J Do JdDi 

Note also that if £, £ H^{D) is such that ^'^^ = V\Di ~ c is constant in Di and 
^'^^^ = Cl-Do harmonic in Dq, then, ^ — cxdi- 

We can decompose uq into the harmonic extension of its constant value in 
Di, c(uo), plus the remainder uq^q € H^[Dq). Thus, we write, 

Uq = Uo.O + c{uo)xDi 

where mo,o G H^(D) is defined by u^l — in D^^' and Uq^q solves the following 
Dirichlet problem, 



/ V4°0-V;^= / fz yzeHliDo). (19) 

JDo JDo 

"o°d ^0 on dDi, 

u'^l ^g on dD. 
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From equation and the observations in Remark [T] we get that 
/ V(uo,o + c(uo)xDi) • VxDi = / /XDi 

J Do Jd 

from which we can obtain 

= JZW^^ ■ ^'"^ 

There is a useful alternative expression for c(uo) in ((20)) that we also use. 
By using the Neumann problem related to mo,o we have that, 

/ Vuo.o • VxDi = / fxDi + / Vuo,o • noXDi 

J Do J Do JdDa 

= / fXDi + / Vuo,o • nol 

J Do JdDi 

and then noting that fxDi = /d„ /xdi + /d, / we get, 

c(mo) = f , (21) 

JdD, ^^D, ■ no 

which reveals that c{uq) balances the fluxes across dDi. To summarize the 
results obtained to this point, we can express uq as follows: 

^ Id fXD, - /g„ Vmq.o • Vxji 
Uo = "0.0 H 7 — 7^ ^ XDi (22) 

= "0,0 + — ^7 XDi (23) 

JdD, ^XDi ■ no 

Given the explicit form of wq, we use it in (jl3p to find u^^^ = from if 

we conclude that Mq^-* and u^^"* satisfy the local Dirichlet problems 

/ Vuo-Vz= / /z VzeH^iDo), 

J Do J Do 

Viii -Vz^ [ fz Vz G Hl{Di), 

Di J Di 

with given boundary data on 9I?o a-nd Equation (jl3p also represents the 

transmission conditions across dDi for the functions u^q^ and u"^\ This is easier 
to see when the forcing / is square integrable. From now on, in order to simplify 
the presentation, we assume that / G L'^{D). If / G L^{D), we have that Mq"^ 
and u^i^ are the only solutions of the problems: 

/ V4°'-Vz=/ fz+ I V^^'-noz Vz G iJ^(L»o) with z = on aZ) 

Jdo Jdo JdDo\dD 
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with u^Q^ = g on dD, and 

/ Vu[^^-Vz^ fz+ VuJ^'-niz ior all z e H^{Di) 

JDi J Di JdDi 



Replacing these last two equations back into (|13p we conclude that 
/ (V4°^ • no + V^^^ • ni)z ioY aW z H^{D), 

JdDi 



which implies 



Vu^i ^ ■ ni = — Vmo"' • no on dDi. 



Using this interface condition we can obtain u\^^ in Di by writing 



ci where / u'{^ = 



and u\^^ solves the Neumann problem 



/ Vu[^^-Vz= fz- v4°^-niz for all z e 

Jdi JDi JdDi 



(24) 

The constant ci will be chosen later. Problem ([24)) satisfies the compatibility 
condition, 

Vm^^^ ■ ni = — I ^"^0°^ ■ "-0 

,(0) 



dDi JdDi 



Vu^ -no- c{uo) / Vxni • no = / /. 

dDi JdD JDi 

Here we use the value of c(mo) is given in (|2T|) . 

Next, we discuss how to compute u^^^ and w^"'' to completely define the 
functions ui G H^{D) and ui £ H^{D). These are presented for general i > 1 
since the construction is independent of i in this range. 

4.1.3 Term corresponding to ?7~* with i > 1: 

For powers of l/t] larger or equal to one there are only two terms in the sum- 
mation that lead to the following system: 



/ Vu,-Vv+ i Vit,+i-Vw = for aU u G i?o(L') 

J Dq J Di 



(25) 



This equation represents both the subdomain problems and the transmission 
conditions across dDi for u^^^ and Following a similar argument to the 

one given above, we conclude that uf^ is harmonic in Do for all i > 1 and that 
is harmonic in Di for i > 2. As before, we have 

Vu|+\ • ni = -Vuf^ ■ no. (26) 
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We note that ul^-* in Di, (e.g., iS^"^ above) is given by the solution of a 
nan 

Uj-^'' = Ci and write 



Neumann problem in Di. To uniquely determine u^'^\ we impose the condition 



W ^ + where / = 0. (27) 

Jdi 



where the appropriate Cj will be determined later 

Given in Di we find 
known Dirichlet data, that is 



Given in Di we find u[^^ in by solving a Dirichlet problem with 



/ Vm|"' • Vz = for all z e iJoi(i:>o) 

^i?o (28) 
ul"-* = itl^-* (= Uj-^-* + Ci) on (9Z?i and Ui = on dD. 

Since c,;, i = 1, . . . , are constants, their harmonic extensions are given by 
CiX-Di J * = li • ■ • ; see Remark [H Then, we conclude that 



Ui=Ui + CiXDi (29) 



where is defined by (|28p replacing c,; by 0. This completes the construction 
of Ui. 

Now we proceed to show how to to find in Di. For this, we use psp 



and ([26| which lead to the following Neumann problem 

[ Vm|+\-Vz = -/" Vwf^-noz for all z G /^^(-Di). (30) 

JDi JdDi 

The compatibility condition for this Neumann problem is satisfied if we 
choose 

lap, - "0 _ Ido - ^^Di /o1^ 

For the second equality see Remark [1] below and Equations ([20|) and ((2T|) . The 
compatibility conditions trivially satisfy 



r7/~(0) , (0)^ 

Vilf ) ■no~cJ Vxg] • no 



= 0. 

where we have used the definitions of Ci given in (1311) . 
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We can choose u. 



(1) 



in Di such that 



'i+l 




+ Q+1 



where 



Ui+i = 0, 



and, as before, 



so we have the compatibihty condition of the Neumann problem to compute 
u\^^2- See the Equation ([30l). 

4.1.4 Summary 

We summarize the Dirichlet-to-Neumann procedure to compute the terms of 
the asymptotic expansion for it^ in (|9|)- p0p . 

1. Compute uq using formulae (|22l) or 

2. Compute w^^'' in Z^i by solving the Neumann problem ()24|) . Compute w^*^^ 
in £>o solving the Dirichlet problem ([25)1 with z = 1. 

3. For 1 = 2,3,... compute in Di by solving the Neumann problem (jSOj. 
Then, compute uf'^ in Dp solving the Dirichlet problem pS]). 

Other cases can be considered. For instance, an expansion for the case where 
we interchange Dq and Di can also be analyzed. In this case the asymptotic 
solution is not constant in the high-conducting part. Multiple inclusions will be 
consider in Section 1431 

4.2 Convergence in H^{D) 

In this section we study the convergence of the expansion (P|)- P^ . For simplicity 
of the presentation we consider the case of one high-conductivity inclusion. The 
converge results will be extended to the multiple high-conductivity inclusions in 
Section [43l We assume that dD and dDi are sufficiently smooth, see (25l |. 

Lemma 2 Let uq in (|22[) . with uo,o defined in (|19p . and ui be defined by (|24p 
and (j28p with i ~ 1. We have that, 



(32) 



\\ui\\hhd,) ^ II/IIh-i(Di) + |l.g||ffi/2(£,) 



(33) 



and 



||wi||h1(Do) ^ \\ui\\m/2(^aD,) \\ul\\H^D^)■ 



(34) 
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Proof. From the definition of uq.o in (fT9)) we liave tliat 

Using p2p we fiave that 

\uo\hHD) = |wo|h1(L>o) < l"0,o|ffi(Do) + |c(wo)||XDilffi(£'o) (35) 
and we observe that 



|c(-Uo)||xDi|ifi(Do) 



IXDilH^Do) 
^ ll/ll-H-i(D) + ||5||_H-i/2(a£))- 



This proves ([5^ . Equation follows from the classical estimate for prob- 
lem (p4)) . Equation (p4)) follows from problem p8| with z = 1 and a trace 
theorem; see 

The following lemma can be obtained using orthogonality relations of Galerkin 
projections. 

Lemma 3 If w E H^{D), w is harmonic in Dq and we define 

w ^ w + c{w)xdi 
where ^ ^ 

c{w) ~ — -T — = — — , 

Jc^D^^XD^ -no J^o I^Xl^i T 

then, w and xDi orthogonal in the operator norm induced by the Dirich- 
let operator, that is, VwVxdi = 0. We also have = + 

c{wf\XD,\Hi{D)' 

\w\m(D) < \w\m{D) and \\w\\hi{d) ^ \\w\\m{D)- 

Here, the hidden constant is the Poincare-Friedrichs inequality constant on D. 

The next lemma bound the norm of the i— th term by the norm of the (i — 1)— th 
in the asymptotic expansion 

Lemma 4 Let Ui defined on Dq by (|28p with q = 0, and w^+i defined on Di 
by ([50)) . For i > 1 we have that 

\\Ut+l\\HHD) ^ \\Ut\\HHD„)- 

Proof. Let i > 1. Consider Ui+i defined by the Dirichlet problem (|28p . From 
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classical estimates of the solution on Dq and the trace theorem on Di, we have 

\\ui+l\\H^Do} ^ ll^»+l|lifi/2(a_Di) ^ \\ui+l\\H^Di)- 

By considering the problem pop we conclude that 

Wc have, form ([3T|) and Lemma |3l we have 
Combining this last three inequalities we have 

\\Ui+l\\H^(D) ^ \\ui+l\\H^(D) ^ IIUjIIhMOo)- 

The constants are independent of i and depend only on the domain geometry 
and configuration, that is, on Di and Dq. In fact, the hidden constants depend 
on the trace theorem and solution estimates in Di and Dq, see (ist . 



Theorem 5 There is a constant C > such that for every rj > C , the ex- 
pansion ([9]) converges (absolutely) in H^{D). The asymptotic limit uq satisfies 
problem (jl4|) and uq can be computed using formula ([22|). 



Proof. From Lemma |4] applied repeatedly i — 1 times, we get that for every 
i > 2 there is a constant C such that 

\\ui\\m(D) < C\\u.i^i\\hi{Do) < C\\ui-i\\m{D) 

< ■ ■■ 

and then 

oo 1 1 1 1 / 

\\2^v u,\\HiiD)< 7^ — 2^1- 

i=2 i=2 ^ ' 

The last expansion converges when rj > C. Using ([32)) and ([33| we conclude 
there is a constant Ci such that that 



1=0 i=0 

Moreover, the asymptotic limit uo satisfies problem p^ . 



Corollary 6 There are positive constants C and Ci such that for every rj > C , 
we have 

I °° f C\ ^ 

\W-Yv''u,\\H^l^D)<Ci{\\f\\H-i(D,) + \\9\\HU2(D)) (-) ' 

1=0 i=/+l \ ^ 

for / > 0. 
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We note that in the case of smooth boundaries dDi , dD and smooth Dirich- 
let data and forcing term, we have and H^~^''{Do) regularity of all 



functions involved for s > 0; see [1^ |25| and references therein. Estimates 



similar to the ones presented in this section will warrant that for 77 sufficiently 
large, the expansion ([9|)-(fT0l) will be absolutely converging in H^~^''{Di) and 
7?^+*(Do) for 77 sufficiently large. A more delicate case is the case with non- 
smooth boundaries. This case and the convergence of the expansion in H^'^'^ [D) 
for some small r > will object of future research. 

4.3 Multiple high-conductivity inclusions 

In this section we consider a coefficient with multiple high-conductivity inclu- 
sions. Let K be defined by 

. ^ r 77, x&D„,, m = l,...,M, 
'^^'^^^ I 1, x&D, = D\ijfUDn.. ^^^^ 

and denote by u^, the solution of (j4]) with zero Dirichlet boundary condition. 
We assume that Di is compactly included in the open set D \ ^JgL-^ ^-^^Df, i.e.. 
An C Z? \ U^£i^„ A, and we define D^:=D\ uf^,^iD„,. 

Expansion (|9|)- (fT0| holds in this case. We first describe the asymptotic 
problem in the next Section 14.3.11 Then we will quickly describe the expansion 
in Section |4 . 3 . 2 1 below . 

4.3.1 The solution of the asymptotic problem 

Define the set of constant functions inside the inclusions, 

Vconst ^ {v E H^{D), such that u|li„, is constant for all m = 1, . . . , A/}. 

By analogy with the case of one high-conductivity inclusion, the asymptotic 
solution for the coefficient (j36]) is uq that is constant in each high-conductivity 
inclusions. Moreover, uq solves the problem, 

Jjy^ Vuo • Vz = fz for aU Z e Vconst, 

(37) 

Uq — g on dD. 

The problem above is elliptic and it has a unique solution. For m = 1, . . . , M, 
define the harmonic characteristic function xd„^ G Hq{D) by 

X_D„. = S„ie on Di for £ = 1, ... , M, 

and, in Dq, XOm is defined as the harmonic extension of its boundary data in 
Do, i.e.. 



/ VxD,„ • Vz = for aU z e H^{Do) 

JDn 



'Do C38-\ 
XDrr, = ^nii on dDi for £ = 1, . . . , M, ^ ^ 

Xd^ ~ on dD. 
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Here, Smi represent the Kronecker delta, which is equal to 1 when m = I and 
otherwise. Remark [1] holds if we replace the one inclusion case xdi with 
the multi- inclusion case xD[ defined in ((55)) . For instance, if u; S H^{D) is 
harmonic in Dq and constant w = c,„ in -D,„, to = 1, ... , A^, then, we can write 

Wc decompose uq into the harmonic extension (to Dq) of a function in Vconst, 
plus a function, uo,Oj with 5 boundary condition on dD and zero boundary 
condition on dDm, to = 1, ... , M. We write, 

A/ 

uo = "0,0 + X! '^™("o)xn„> (39) 

m— 1 

where uo,o G H^(D) with uo,o = in Z3„i, to = 1, . . . , M, and wo,o solves the 
following problem in Dq, 

/d Vuo.o • Vz = /z for all z e iJ(}(Do), 

(40) 

1*0.0 = on dDm, to = 1, ... , A/, and uo,o ~ 5 on 

Equation (|39p is the analogous to Equation (|22p . Now we show how to compute 
the constants Ci(uo) using the same procedure as before. From (jl4p . we have 

A/ 

for ^ = 1 M, 



/ V(uo,o + V c™(uo)xn™) ■ VxDf = / fXDc, 



which is equivalent to the M x M linear system, 

Ag^omX = B (41) 
where A = [a^], and B = (&i, . . . , 6^/) £ R^^ are defined by 

a^J= I VXD. VXD, = / VXD. VXD,, (42) 

JD J Da 

bj = / fXDj - / Vwo.oVxd, 

J D J Do 

and X = (ci(uo), . . . , C7\/(uo)) G K^^. We conclude that 

X = A-l,^B. (43) 
Wc note that using ([T7)) for XDi wc have that 

= / VxDiVxD, = / ^XDj-rii^ / VxD, -nj- (44) 

JD JSDi JSD, 

Note that X]m=i Cm(uo)xDm is the solution of a Galerkin projection in the 
space Span{xD™}m=i- The forcing term for this problem is / and there is 
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Neumann boundary data on dDm coming form S/uofi-no- Matrix Ageom encodes 
the geometry information concerning the distribution of the inclusions inside the 
domain while it is independent of the contrast 77. Note that in general Ag^o^ 
can be a large dense matrix. Because xOj decay, one can approximate the system 
by a sparser system (e.g., see Moreover, we can use concepts similar to 

multiscale finite element methods and seek smaller dimensional approximations 
for this large system. 

4.3.2 Expansion 

Now we describe how to compute the individual terms of the asymptotic expan- 
sion (P))- ([T(I)) for the case of multiple high-conductivity inclusions. 

• The function uq solves (|37l). 

• The restriction of ui to the subdomain Dm, Ui^\ can be written 

wj™) + ci,„ where / Sf^™) = 0, 

and Wj™^ satisfies the Neumann problem, 

/ Vu["'^-\/z= I fz- I Vu^Q^-n^z for all z e i7i(L»™), (45) 

for TO = 1, ... , M . The constants ci^,„, m = 1, . . . , M, will be chosen later. 

• For i = 1,2,..., we have that given uf^^ in Dm, to = 1, ... , M, we can 
find uf*^ in Dq by solving the Dirichlct problem 

( Vuf ^ • Vz = for aU z G HI{Do) 

J Do 

(46) 

= uj™^ (= "i™^ + Ci,m) on SAn, TO = 1,...,M, and 
uf^ on dD. 

Since Ci^m are constants, the corresponding harmonic extension is given 
by Hm c-i,mXD( ■ Then, we conclude that 

= Mi + ^ Ci,™X-D™ (47) 

m 

where uf''^ is defined by (|46|) replacing all the constants Ci,m by 0. 
The u[™l in D„i satisfy the following Neumann problem 

/ Vu|"i^ • Vz = - / Vuf ^ ■ noz for aU z e H\Di). (48) 
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For the compatibility condition we need that for £ = 1, . . . , M, 



= / Vu|+\-nf 

JdD, 



dDe 

M 



f M 

From (|42p and we conclude that "K^ = (ci^i, . . . , q^a/) is the solution 
of 

where 

[/=(-/ VSff^.no,...,-/ Sff^.no). 
or (using Remark [1]), 



U={- 



4.3.3 Convergence in H^{D) 

We first prove the result analogous to Lemma |3l 

Lemma 7 Letw G H^{D) he harmonic in Dq and define w = 'S'+X)m=i CmXn„ 
where Y ^ (ci, . . . , Cj\f ) is the solution of the M dimensional linear system 

AgeomY = —W 

with W = (/^^^ Vu; VxDi , • • • , /o„ V«;VxD„ ) • Then, 

\\w\\m(D) < \\w\\m(D) 
where the hidden constant is the Poincare-Friedrichs inequality constant of D. 



Proof. Note that CmX-D™ is the Galerkin projection of w into the space 
span{x_Di}f£i- Then, as usual in Finite Element analysis of Galerkin formula- 
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tions, we have 



» M 

J Dn „,_i 



E 



» M 

/ V«;V(Ec™XD„) 

•'^o rn=l 
M 



< \Mm{Do) 



E 



m— 1 



and then 



X]m=i '-'TiX-Dm ^ \'^\h^(Do)- Using a Poincare-Friedrichs hiequality 



we can write, 



\\w\\miD) < \\w\\h^D) + 



M 

E "^'"XDn 



HHD) 



Combining Lemma [7] with rcsuUs analogous to Lemmas [5] and U] we get 
convergence for the expansion (|9))- (|T0)) . 



Theorem 8 Consider the problem (|4]) with coefficient (j36p . T/ie correspond- 
ing expansion ([9 |) -(|10 p converges absolutely in H^{D) for rj sufficiently large. 
Moreover, there exist positive constants C and Ci such that for every rj > C , 
we have 



I 

E 

■i=0 



/ C 

V^'u^\\hi(d) <Gi{\\f\\H-i(D^) + \\g\\mn(D)) E (7) ' 

i=I+l ^ ^' • 



for / > 0. 



5 The case of low-conductivity inclusions 

In this section we derive and analyze expansions for the case of low-conductivity 
inclusions. As before, we present the case of one single inclusion first (see 
Section [Q]) and analyze the general case in Section [5731 



5.1 Expansion derivation: one low-conductivity inclusion 

Let K be defined by 



K,{x) 



e, X & Di, 

1, X (E Do = D\Di, 



(50) 
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and denote by Ue the solution of We assume that Di is compactly included 
in D {Di C D). Since Ue is solution of Q with the coefficient ([50)) we have 

/ Vu, -Vw + e/ Vu,-Vv= / fv for sdl v e H^{D). (51) 

Jdo J Di J D 

We try to determine C Hq{D) such that, 

oo 

= e^^u^i +uo + eui + e^U2 + • • ■ = ^ e'uj, (52) 

and such that they satisfy the following Dirichlet boundary conditions, 

"0 = .9 on dD and ui = on dD for i = — 1, and t > 1. (53) 

Observe that when w_i ^ 0, then, does not converge when e 0. 

If we substitute ([5^ into (|ST|) we obtain that for all v G Hq{D) we have, 

/ Vwo-Vw + VeM / Vu, ■ Vu + / Vuj_i-Vw) = / /w. 

JDo i^o ^-^Da JDi J Jd 

Now we equate powers of e and analyze all the resulting subdommain equations. 

Term corresponding to : 

We obtain the equation 



/ Vm_i • Vw = for all w e -ffoP)- 



(54) 



Since we assumed w_i = on dD, we conclude that V?i_i = in Dq and then 
= in Da. 

Term corresponding to = 1 : 

We get the equation 

/ Vuo-Vw+ / Vu_i-Vu=/ /w for all w e iJo^(L>). (55) 

Jdo J Di Jd 

Since u"^\ = in Dq, we conclude that ul\ satisfies the following Dirichlet 
problem in Z?i, 

/^^ Vu^^l ■ Vz = /^^ fz for all z e Hl{Di) 



(56) 

mL^] = on az?!. 

Now we compute Uq"' in Do- As before, from ([55]). 



Vuq"' • no = — VwL^J • rii on 
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Then we can obtain u'j^^ in Dq by solving the following problem 
Xdo "^"o"^ ■ '^^ /do -/"^ " /9D1 ^"-^1 • "1-^ ^ ^^^i^o) with z = on dD, 



4°^ = ,g on dD C aZJo. 



(57) 



Term corresponding to with i > 1: 

We get the equation 

A) "'-Di 

which implies that u'"'^' is harmonic in Di for all i > and that w^*^' is harmonic 
in Do for i > 1. Also, 



V7 (0) V7 (1) 

Given u\^^ in Z?o (c-g-, wq in Dq above) we can find in Di by solving 
the Dirichlet problem with the known Dirichlet data, 

Vu^i\ • Vz = for aU z e H^{Di) 

(58) 

(1) (0) cm 

To find uf^'' in Dq we solve the problem 
/ Vwf ^ ■Vz = - I Vw,|i\ • niz for all z S ^^(Do) with z = on ^D, 

Jdo JdDo 

(59) 



= on dD. 



5.2 Convergence in H^{D) 

In this section we study the convergence of the expansion ([5^ -([5 ^ . The fol- 
lowing lemma is obtained using classical estimates and trace theorems in the 
involved subdomains. 

Lemma 9 Let it_i vanish in Dq and he defined using problem (j56p in Di. We 
have 

1^4-1 ||_H"i(-Di) ^ 

and 

\\uo\\hhd„) di WIWh-^Do) + \\U-i\\hHDi) + hWH^^dD)- 
Moreover, if we consider problems (j58p and (|59p . we have that for i > 1, 

\\Ut\\H^Di) d ll^illifMA))' \\ui\\H^D„) ^ \\ui-l\\H^Di) 

\\ui\\H^{D) di ||wi-l||_ffi(-Di)- 
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The convergence of the expansion foUows. 

Theorem 10 There is a constant C > such that for every e < 1/C, the 
expansion (j52p converges (absolutely) in H^{D). 



Proof. There is a constant C such that, for every i > 1 we have 

\\ui\\m{D) < C\\ui^i\\Hi{Di) < C\\ui^i\\m(D) (60) 

< ... (61) 

< C'WuoWhhd,) (62) 

and then 

\\}_^eu.,\\Hi{D) < ^ 2^ (Ce) . 

1=1 i=l 

The last series converges when e < 1/C. Using the bound for uo,ui and m_i 
we obtain that there is a constant Ci such that 



E^'^^II^Mi^) < Ci{\\f\\H-HD,) + II/IIh-(d,) + \\9\\mmaD))J2^^'y 

i=0 i=l 



Corollary 11 There are positive constants Ci and C such that for every e < 
1/C, we have 



I 

\\u ~ u 

i=l 



1 -^e'uiWH^D) < 



CimH--iD,,) + \\f\\H-HD^) + U\mn(OD)) E (^^)' 



for / > 0. 



Using the asymptotic expansion (j52p - (j53p we can write an asymptotic prob- 
lem for e — !■ 0. 

Corollary 12 If f — Q in Di, we can write 

We = uo + eul + e^U2 + ■ ■ ■ 

where u^^^ = uqIdq satisfy the following problem with Dirichlet data on dD and 
zero Neumann data on dDi . 

/^^ V4°^ ■ Vz = fz Vz e HI{Dq) with z = on dD, 
Ug^' = g on dD C dDg. 
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Additionally, we can find Uq = Uo\di by extending harmonically to Di the 
known Dirichlet data on dD, that is, 

[ Vu^\-\/z = for all z e H^{Di) (63) 

with 

(1) (0) on 

Uq ~ Uq on oL>i. 
The series converges absolutely in H^{D) for e sufficiently small. 




Figure 2: Example of disconnected background. Here Dq ~ AU C and Di = B 
is a low-conductivity inclusion. Sec Remark 1131 

When Dq is not connected, the following observation can be made. 

Remark 13 In the case of Dq being disconnected, we have that (|54|) implies 
that w_i is constant in each connected component of Dq and it vanishes only 
in the connected components whose boundary intersects dD. The function u_i 
will be constant in the other interior connected components. This is similar to 
the case of high- conductivity inclusions. The function U-i will be zero only if 
the forcing term vanishes in these interior components also; see Equations pOp 
and (|49|). For instance, consider the case illustrated in Figure [H where the 
background domain are Dq ~ AU C and the inclusions is given by Di = B. In 
this case it is easy to see that u_i satisfies a problem similar to problem (jl4p 
in B U C. Then, u_i will vanish only if the forcing term vanishes in B U C. In 
this case, a result similar to Corollary [T^ above can be stated. 



5.3 Multiple low-conductivity inclusions 

Let K be defined by 

J e a: e An, m = l,...,M, 
"("^ = 1 1, ^eDQ^D\U^^^,D^, 

and denote by Ue the solution of ([4]) with coefficient (j64| . We assume that Di is 
compactly included in the open set D\\jf^i^^^j^De, i.e., Dm C D\'jfLi g^j^De, 
and we define Dq:= D\ uf^^Jj^- 
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Expansion extends easily to this case of multiple low-conductivity 

inclusions, that is, 

• Wc have u^^J = on Dq- Also, that each, satisfies the following 
Dirichlet problem in Dg, 

J^^ ■ Vz = J^^ fz for aU z e Hi {D,) 

(65) 

uL^J =0 on dDi. 

• We can obtain Uq*^' in Dq by solving the problem 

M 



[ Vuf-Vz^f fz-f^f Vu^'\-niz ^zeH^iDo),z\9D = 0, 

J Do J Do JdDe 

Mo°^ = g on dD. 



(66) 



• Finally, we have that u\^^ is harmonic in Di for all i > and that uf^ is 
harmonic in Dq for i > 1. 

Given uf^^ in Dq, we can find u^j^2i in Dg by solving the following Dirichlet 
problem with the known Dirichlet data, 

/ Vuf \ • Vz = for all z £ HUDA 
Jd, (67) 

ufli — uf\ on dDi. 
To find uf''^ in Dq we solve the problem 

/ Vu'-p . Vz = - V / Vuf2i ■ UiZ Vz e HI{Dq) with z\qd = 0, 

J Do «_n JdDi 



I dDi 

Ui=Q on dD C dDo. 



(68) 



The convergence of the expansion (|52p - ([55)) is similar to the case of one low- 
conductivity inclusion. In particular Corollary [TT] holds in this case. 



6 An example with low- and high- conductivity 
inclusions 

In this section we show an example with a high- and a low-conductivity inclusion. 
The procedures were introduced in detail in Sections |4] and [5] We show only 
how to write the subdomains problems for the leading terms of the expansion. 
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Consider k to be defined by 

{77, X e Di, 

xeD2^ _ (69) 
1, X e Do = D\{DiUD2). 

As before we write u,, = 7?m_i + ua + + ;^M2 H = Si^-i 'H'^'^ii witli 

Wo = 5 on and = on for i 7^ 0. We need tliat 

/ Vu-i-Vw + r/ / Vu-i • Vi' + / Vmq • Vv ) + 
^ri'([ Vu,-Vv+ f Vu,+i-Vv+ f Vu^-i-Vv]^f fv. 

i=0 ^•'Do JDi Jd2 / Jd 

Now we equate powers and analyze the subdomain equations. We assume tliat 
Di is connected and compactly included in D {Di C D), i = 1,2. We also 
assume that the distance between Di and D2 is strictly positive, then u^^l = 0, 
u^^l = and u*^] solves the following Dirichlct problem in _Di, 

/^^ ^u'-ll ■ Vz = /^^ fz for all z e H^{D,) 

(70) 

uL^J = on dDi. 
This defines the function 7i_i. To write a problem for uq, let 
Vconst ~ {v G {D\D2), such that w = on dD and u^^-* = uId^ is constant }. 
We have that 

/dh • Vz = JdJz- Jqd^ VwL^J • for aU Z G Vconst, 

(71) 

Uq = g on 91?. 

The problem above can be analyzed using the harmonic characteristic func- 
tion xdi G Hq{D) defined in ([55)1 with M — 2. The solution of the asymptotic 
problem above gives Uq°'' and the constant function u^^K As before, the constant 

Uq^^ can be determined explicitly using an expression similar to (|2ip. To com- 

(2) 

plete the definition of uq we observe that Mq ' satisfies the following Dirichlct 
problem with the known Dirichlct data, 

Vu'q '^ • Vz = for all z G H^{Di) 

(72) 

(2) (0) an 

Uq — Uq on oLli. 
The functions Ui, i ~ 1, . . . , can be determined form the equation 

/ Wui-\7v+ Vu,+i-Vv+ Wui^i-Vv = for all V e Hq{D). 

J Do J Di J D2 
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This procedure is similar to the ones developed before and presented in detail in 
Sections |4] and [5] As before, Ui is harmonic in each region. Its restriction to sub- 
regions can be determined by solving subdomain problems involving Dirichlet, 
Neumann or mixed boundary conditions on the inclusions boundaries. 



7 Conclusions and comments 

We use asymptotic expansions to study high-contrast problems. We derive and 
analyze asymptotic power series for high-contrast elliptic problems. We mostly 
consider the case of binary media with interior isolated inclusions. High- or 
low-conductivity inclusion configurations are considered. The coefficients in the 
expansions are determined sequentially by a Dirichlet-to-Neumann procedure. 
In the case of high-conductivity inclusions, the Neumann problem needs to 
satisfy a compatibility of fluxes. This flux-compatibility condition is obtained 
using an auxiliary finite dimensional projection problem. The related finite 
dimensional space is spanned by harmonic extension of characteristic functions 
of each subdomain. 

The asymptotic limits when the contrast increases to infinity are recovered 
and analyzed; see Theorems [5] and [8] and Corollary [121 The convergence of 
the expansions in H^{D) is obtained provided that the contrast is larger than 
a constant C that depends on the background domains and the domains rep- 
resenting the inclusions. The convergence rate is algebraic. We consider the 
case of isolated interior inclusions which can be high and low conductivities. 
Other more complex configurations can be analyzed following a similar proce- 
dure. The analysis covers the cases with high- and low-conductivity inclusions. 
See Figure [3] for schematic representations of two dimensional configurations. If 
the low-conductivity value is e and the high-conductivity value is 77 = 1/e, then, 
an expansion similar to ((52)) -([53 )) can be used for all the examples in Figure [3l 
For general values of 77(e) = e~P, an expansion similar to can also be 

derived where the coefhcients in front of spatial terms will scale as e'+™P, where 
/ > — 1 and ?Ti > are integers. 

More general coefficients can also be studied. Similar expansions for other 
problems related with flows in high-contrast multiscale media can be obtained, 
e.g., models like heat conduction, wave propagation, Darcy or Brinkman flow, 
and elasticity problems. Efficient solution techniques for solving the system 
of linear equations (PT|) will be a subject of future research, in particular, lo- 
calization procedures for the harmonic characteristic functions will be studied. 
Questions concerning the convergence of the series in stronger norms as well as 
computing quantities of interest will be studied in the future. Reduced contrast 



approximation and related multiscale methods as in [12|, [13| will be the subject 
of future studies. 
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Figure 3: Four different configurations witli liigh-conductivity 1/e (red) and 
low-conductivity e (blue) regions. Expansion (|52p can be obtained in these 
cases. 
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